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Abstract 

This paper develops further approximate methods for obtaining the dipole 
matrix elements and corresponding transition and decay rates of the high- 
n, high-/ gravitational eigenstates. These methods include (1) investigation 
of the polar spreads of the angular components of the high-«, high-/ 
eigenstates and the effects these have on the limiting values of the angular 
components of the dipole matrix elements in the case of large / and m and 
(2) investigation of the rapid cut off and limited width of the low-/», high-n 
radial eigenfunctions, and the development of an equation to determine the 
width, position and oscillatory behaviour of those eigenfunctions in cases 
of arbitrarily large values of n, I and p. The methods have wider 
applicability than dipole transition rate estimates and may be also used to 
determine limits on the rates for more general interactions. Combining the 
methods enables the establishment of upper limits to the total dipole decay 
rates of many high-«, low-p states on the state diagram to be determined, 
even those that have many chaimels available for decay. The results 
continue to support the hypothetical existence of a specialized set of high- 
n, low-/? gravitational eigenfunctions that are invisible and stable, both 
with respect to electromagnetic decay and gravitational collapse, making 
them excellent dark matter candidates. 

PACS numbers: 03.65.Ge, 03.67.Lx, 03.65.Db, 04.60.-m, 95.35.+d, 
04.90.+e 



1. Introduction 

In a previous paper [1], a study of the properties of gravitational eigenstates was 
begun, the motivation for which originated from earlier work [2] which 
suggested that certain members of the quantum eigenstate ensemble of deep 
gravitational wells might have properties that may make them relevant for 
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understanding the nature and origin of some astroptiysical phenomena, such as 
Dark Matter. The motivation for further theoretical study of gravitational 
eigenstates is made even more compelling by the recent experimental 
verification of their existence [3, 4]. 

In particular, paper [1] developed general, non-integral formulae for calculating 
the transition rates of state-to-state dipole-allowed transitions vwthin the 
eigenspectra of gravitational potential wells. Sometimes these non-integral 
formulae continue to have limitations caused by a large, non-truncatable number 
of summation terms, and/or excessively large numbers, and other more 
approximate methods need to be employed. The purpose of this paper therefore 
is to develop further practical methods for these cases which are of use in 
approximating matrix elements for both the spontaneous dipole decay of the 
high-quantum-number gravitational eigenstates, and potentially also of more 
general interactions. In particular we examine ways of determining the 
approximate overall position, size and internal structure of the individual 
eigenstates, in particular their radial components. 

The position, shape and structure of the eigenstates turn out to be very useful in 
determining certain dipole transition overlap integrals and can aid in setting 
maximum limits on the absolute value of others. They are also important in 
investigating limits on higher order (multi-pole) decays and overlap integrals for 
other more general interactions, with both electromagnetic radiation and 
particles. (For example, it is clear that if two states share no common region of 
space then no interaction Hamiltonian could induce transitions between them.) 
Structural features of the wavefunction components of an overlap integral can 
also strongly influence its value even when the states do share a significant 
spatial overlap region. Indeed there are particular individual eigenstates whose 
structure and position are such that the only significant transitions possible are 
those that have a negligible effect on the perturbing entity. 

Important structural features of eigenflinctions include variation of the average 
absolute value of the wavefunction across its spatial range, and the spatial 
oscillation frequency (SOF) of the zeros or nodes of the eigenfimction 
components. We wish to establish how these features vary with the location of 
the state on the quantum state diagram. This diagram is reproduced here from [1] 
as figure 1 , and uses the quantum parameter p[=n-l) as defined in that paper. 



2 



Gravitational eigenstates in weak gravity II: 



l=n-2 l=n-4 



increasing / value 



/ = 




increasing 
n value 



Figure 1. Schematic representing the high n, I , m -valued stationary 
states, drawn to emphasize the parameter p = n — l . Each solid circle 
on the diagram represents 2/ + 1 z -projection substates. 

More specifically, for a given eigenstate on figure 1 , we wish to know, how far 
away on figure 1 a second state must lie (in terms of changes in the quantum 
parameters Ah, A/ and/or Aj?) so that differences in the structural features of 
the two states will ensure a negligible probability interaction. This will depend 
on the location of the initial state and on the direction of the second state. 
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In section 2 we consider ttie effect that the size and shape of the angular 
component wavefunctions have on the value overlap integral, particularly in the 
limiting case of large quantum numbers, and derive limiting values for the 
angular component of the dipole matrix element. Section 3 deals with the 
corresponding radial component of the dipole matrix element. We firstly look at 
exact results for very low-p that demonstrate the rapid reduction in the radial 
overlap component for high-« as the difference between the quantum state initial 
and final p values Ap and corresponding principle quantum number change An 
increases. We then examine a procedure developed to determine the approximate 
positions of the zeros of the radial eigenstate wavefunctions, which enables the 
width, spatial oscillation frequency and average position of any radial 
eigenflinction to be obtained for arbitrary values of n and /. Combining the 
results from section 2 and 3 then enables decay rates for various significant 
quantum states to be obtained and this is done in section 4. Section 4 also deals 
with calculating the decay rates of other constant-/) transitions, showing that for 
low-/», the rates are similar to those for the p = \ states. In section 4.3 transitions 

involving A« = 2,3,4..., but those that do not necessarily end on p-\ are 

examined, and procedures and formulae for dealing with the \ow-p states with 
multiple decay channels is developed. This section then uses some of the 
methods developed to obtain trends in decay rates as one moves vertically and 
horizontally away from the p-l diagonal of figure 1 , demonstrating how the 
long hfetimes of states decrease as one moves away from the low-/* diagonals. 
We conclude with a summary and discussion of the salient points of the paper in 
section 5. 

2. The effect of the angular eigenfunctioii components on 

We consider separately the spherical harmonic (angular) and associated 
Laguerre (radial) eigenfunction components. As in [1] the total dipole decay rate 
is given by 



where as in [1] e is the electronic charge, ^^e electrical permittivity of free 



corresponding to the transition z to / , // the reduced mass, and the other 




(1) 




symbols have their usual meanings. is given again as in [1] by 
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^ifA\o\l\,^^nfffYirnferRniAmir'^^^\0)^om dct>d0dr 

\ (2) 

where 

h<j,x^\l\l'' Yi}^^f \^iSm{0)cos{^)sm{0)d^0 , 

^o<l>y = Jo Jo ^VM \mi sin(^)sin(^^)sin(^)J«^f/6' and 

V = \l\T ^VM \mi cos(0)sm(0)d^0 
with the R„j ii being related to the associated Laguerre polynomials and 7/,^, the 
spherical harmonics. 

2.1. The effect of shape and spread of the probability density distribution on IT^y 



The azimuthal {<j)) dependence of all spherical harmonic eigenfunctions 

components is of the form e*' " so that the probability density in the azimuthal 
direction is uniform for all states, and hence there is always functional overlap in 
this coordinate direction. The spherical harmonic polar {0) dependence varies 
with the values of both / and m. A zero angular momentum z-projection value 
(i.e. m = 0) corresponds to the highest degree of polar angular spread for a given 
/ value, with the lobes of probability density equally spaced around the polar 
direction (figure 2, curve (/)). A state with a high z-component of angular 
momentum (m close to / ) however corresponds to a polar distribution that is 
more flattened and disk-like, but with the probability density lobes still 
approximately equally spaced and centred around 0-7i: (figure 2(a)). Clearly 
for any two spherical harmonics, there will always be some degree of overlap in 
the polar direction as well and the effect that this has on I^g^ , I^g^ and g,^ in 

the limit of large / needs to be determined. It also means that two eigenstates 
can be spatially disjoint only if the radial components of their wavefunctions are 
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physically separated in space. We will return to the radial components after 
examining the limiting values of the angular part of the dipole transition overlap 
integral. 
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Figure 2. Polar angular position (0 to tt) of the zeros of the 
probability density functions of the spherical harmonics as a 
function of zero number for / = 200 and various m values, (a) m = 
180; (b) m = 160; (c) m = 120; (c/) m = 80; (e) m = 40; (/) m = 0; 



2.2. Limiting values of the angular components of the overlap integral 

Unlike , which can become arbitrarily large because of the presence of the 

extra factor of r, the integrals Iq^^ , Ig^y and Iq^^ are limited in size even when 

there is considerable overlap. Since we are examining states where the values of 
n are large but p is small, we will generally be concerned with determining how 
the size of these integrals behaves when the initial value of / is large. The initial 
values of m are arbitrary however and each of the points on the low-p diagonals 
of figure 1 represents all allowed values of m from -/ to +/. There are essentially 
two limiting regimes. Either m remains small as /^oo or both / and m 
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with I ~ m. Table 1 shows the limiting behaviour of angular components Ig^^ , 
Ig^y and Iq^2 , (a) when m is small compared to / and (l-m)^oo, and (b) 

when \m\ is near / and /,m— >oo. Dipole transitions are accompanied by a 
angular momentum change A/ = ±1 and a z-projection change of Aw = ±1 . 

Table 1. Limiting values for the angular component of the dipole matrix 
elements for electromagnetic decay: (a) lim(/ - m -> oo) , when m«l and (b) 

lim(/, m oo) when m-l . 
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2.3. The effect of limiting values of the angular components of the overlap 
integral on eigenstructure shape 

The values in table 1 can be used to predict trends that might occur when upward 
and downward transitions take place (for example due to stimulated emission 
and absorption that could occur in instances where photon densities were 
sufficiently high over the whole of the eigenstate's volume). For states whose z- 
projection angular momentum value m is close to zero (table 1(a)), both upward 
and downward changes in m are equally likely. If sufficient radiation density 
were present to result in significant numbers of stimulated transitions then the 
effect of this would be to diffuse the state population evenly around m = 0. 
However for those states that have an absolute value of m close to /, transitions 
take place in a manner better described by table 1(b). This shows that for 
transitions where / is increased, states with positive m have their m value further 
increased and states with negative m values have their m values further 
decreased to more negative values. Thus transitions to larger / favour a net 
migration towards states with an increased z-component of angular momentum, 
(more 'oblate'). Conversely transitions to lower / are accompanied by a net 
migration toward states of lower z-component of angular momentum, that is, 
more highly spherical states. 

In the atomic case there is no significant correlation of change of n value with 
change of /. However as will be seen later, in high-n states crossing p diagonals 

in a transition has a large effect on the value of , and increases and decreases 
in / are correlated with respective increases and decreases in n. (That is for 
example, a A/ of +1 is much more likely to be accompanied by a An of +1 
rather than a A/7 of —1 , for high n state.) Then, since state size is related to the n 
value, an eigenstructure containing many particles might be expected to exhibit 
differences in shape that reflect the degree of stimulated emission and absorption 
that has taken place in its history. A randomly populated eigenstructure ought to 
be spherical, since the sum over all m states of any fixed total angular 
momentum value / is spherically symmetric. If however an eigenstructure was 
subject to significant stimulated transfers in its past it might have an inner 
structure that is more spherical but an outer structure that could be somewhat 
oblate or flattened. 

3. Radial component 
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We now return to the radial components of the eigenfunctions and the radial 
dipole overlap integral component Ij^ = r R„ijidr which is of most 

importance in determining the radiative lifetime of any high quantum number 

state. We recall the significance of the p value on the state diagram of figure 1, 
that is that it reflects the spatial oscillatory behaviour of the eigenfunction (p= I 
corresponding to a radial eigenfunction with one peak, no zeros; p = 2, two- 
peaked, one -zero; etc.) We are interested in the "beginning" and "ending" 
position of the eigenfunctions as well as their total radial spatial extent and their 
spatial oscillation frequencies because these things have the most profound 
effect on determining the value of the overlap integrals. The radial wavefimction 
components have the greatest radial extent when / = 0, from low r up to 

r~2n Zjq where is a scale factor [2]. 

3.1. Exact calculation of Ij^ and eigenfunction position, size and spread, for 
low-p eigenstates, and the demonstration of rapid decay 



When the / value is very close to n, the Laguerre polynomials contain only a 
small number of terms (determined by the value of p = n-l), and the physical 
structure and overlap integrals can be examined exactly. We define the width of 
a p = \ eigenfiinction as the full width at zero concavity (FWZC). The peak 
centre and FWZC of a p = \ (Rf,,i=„-i) eigenfunction are given by 

n{n-V)bQ ~ n^b^ and In^n-lb^ ~ In^^'^b^ respectively. This means that two p 

= 1 radial eigenfunctions need differ in n values hy ~4n before their radial 
position is significantly different compared to their width. For larger separations 
the amplitude of either function at the position of the peak of the other drops 
very rapidly. It may be shown that given two p = \ eigenfunctions whose peaks 
are separated by w FWZCs, that the ratio of the amplitude of one wavefunction 
at the position of the peak of the other, to the amplitude of the peak of the first is 
given by 



exp[-2w^«;-l] 



exp 



1 + 



2w 



\n,-l 



(3) 



and that the radial dipole overlap integral has an approximate value of 



\5/2 



(4) 



where n, and nf 



are the (» 1) initial and final quantum levels. 
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The dramatic behaviour of (3) and (4) for large n is demonstrated in the 
following example. For a central mass of 10^" kg with n = 10^^ (binding energy 

~2 eV at a radius of ~ 6x 10^' m for a proton mass), the p-l state has only one 

possible downward dipole-allowed transition: Ap-0, An = l. The initial and 

final states exhibit virtually 100% overlap but an angular emission frequency of 

coJuJ 

less than 0.1 //Hz, so that the dipole decay rate is negligible (A^f- — ^ ). 

'is^Kftc 

Multi-pole electromagnetic transitions and other interactions can however allow 
An > 1 transitions so the degree of separation of two = states for arbitrary 
An becomes an important consideration for determining the maximum value of 
an overlap integral in these cases. Taking the dipole decay as an example 
however, as A« increases the overlap gradually decreases and so does /^j . For 
w=l, (3) gives the amplitude ratio as 0.4, but for >v=10 (corresponding to a 
transition with An~3xl0^^) it is 1.4 xlO"^'' and for w=100 it is ~lxl0~^^^^. 
The corresponding /^j values are ~ 3 x 10" for w=l, ~lxlO~^" for w=10 and 

~lxlO~^^^^ for w=100. The overwhelming decrease in the corresponding 
overlap integral when the two wavefunctions are spatially distinct far 

outweighs increases in A^j - — — that might result from increases in other 

parameters due to larger An , such as «y,y in this case (for w= 10 , coif increases 
to ~ 50kHz). This clearly demonstrates the effective extinction of interaction 
between two p-\ states as w increases. 

For larger n the effect is even more dramatic: for a lO"*^ kg central mass 
potential, a transition that takes place from a n = lO^'', p = \ eigenstate, requires 

a quantum level change of An ~ 3x10 to produce a 50kHz photon. However 
the degree of overlap for this value of An is such that the radial component of 

the dipole matrix element Ij^ = J^nf,if ^ ^niji'^^ ^ value of only 

~ 10"^-^''^*' . A similar analytical procedure can be undertaken for p = 2, 3, 4... 

with similar results except that, because the states have a progressively larger 
spatial extent, the required value of An to ensure distinct separation of states in 
space is now be larger. This concept is illusfrated diagrammatically in figure 3. 
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decreasing / 
► 




Figure 3. Schematic diagram showing changes in position of high-n, 
low-/?, radial component eigenflinctions with n, and the significance of p 
in determining the change in n required for radial distinct radial 
seoaration of the eieenfunctions. 



3.2. Position, width and zeros of the high-n radial eigenfunctions 

The importance of the width, position and SOF of any given eigenfunction 
means that it is useful to have a reliable way to determine these structural 
parameters. We now develop an equation for determining the position, extent 
and zeros of the radial eigenfunctions. Figure 4 shows the typical oscillatory 
form of the radial function R„ j r and diagrammatically shows the way that the 

terms "first zero spacing", last zero spacing" and "estimated polynomial width" 
are defined. The SOF is defined as the reciprocal of the zero spacing at any 
position. 
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first zero spacing 




last zero spacing 




estimated polynomial width 



Figure 4, Schematic diagram showing the typical oscillatory 
behaviour of the function i?„,(r) r, the first and last zero spacing 

and the criteria for estimating the Laguerre polynomial width. 
R„i r , written in terms of p, is 



f 2 



\p-\)\{2n^-p)\ 



In, 



exp 



r X 



(-1)' 



y 

^^^[p-k,-\)\{ln,-lp^k,^\)\k^\ 



V 



(5) 



A little algebra shows that this function, and also the functions R^ i and 

\2 



[r^ I , have the same zeros as the polynomial: 

tQ{p-\-k)\{ln-lp^k^\)\k\ 



(6) 



where y-lrj^nh^^ and again p-n — l. Once the properties of ^(7) are 
determined, it is then a straightforward matter to obtain the properties of i?„ / , 

\2 



^nj r , (Rnj r) etc. 
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Figure 5, Square roots of the 1"', 2'"^, 3'^, ... f\ ... zeros of the function 

(^i)*""^' (rY (p - 1)'(2« - n)' 
^(r) = S ^ — — as a function of the zero number. 

Curve set (a): j7 = 20, « = {21,...30}; curve set (Z)): p = 15, « = {41,...50}; 
curve set (c): js-^lO, « = {61,...70} ; curve set (J): p = 5, n = {Sl,...S9} . 



We wish to obtain from (6) estimates of the spatial extent and position of ^{y) , 

the spacing (spatial periods) of its zeros, and the height of its peaks. The process 
begins with an empirical observation shown in figure 5 that the set of square 
roots of the zeros of Y{r) , plotted as a function of zero number, has a consistent 

characteristic form that is approximately linear. The degree to which linearity is 
approached is independent of the value of or « separately but is better when 
p«n and n»\ concurrently (curve set {d), figure 5). Y{y) is a polynomial 
of degree p-\ and may be written in terms of its p-\ zeros. 
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n^Yi ,r3^---rp-\ as '^[r)^{r-r\){r-r2){r-n)--{r-rp-i)- 

Mathematically, if the square roots of the zeros of T(;') are linearly related, 
then the z* zero Y^ may be written in terms of the first zero y-^ as 



(7) 



where ^ is a 'spacing' parameter that is constant for a given value of n and p. 

p-\ 

We let the sum of the zeros be « = and the siun of the products of the 

p-\ f p-\ p-l p-l '\ 

zeros taken two at a time he fi= ^ -~ Yi Yj ~ S 

Substituting -J/i = -Jxi +{^~^)^ and -Jyi = -Jyi into these 
summation expressions for a and P and simplifying gives: 



« = ^(p-l)((2/ -7/^ + 6)^2 +6(/p-2)^7^+6n) 

(20/ -156/ +427p^ -477;? + 180)j^ + 
120(/-6/+ll/7-6)(J^^ + 
60(5/-20/j + 18)<5Vi + 

From simple polynomial theory -a is the coefficient of y''^^ and P the 
coefficient of 7^ in equation (6). Using equation (6) gives 



360^ ' 



(8) 



(9) 



a = {2n-p){p-\) 
P = {2n-p){2n-p-\){p-\){p-2)l2 
Eliminating a and /3 using equations (8 - 11) gives: 



(10) 
(11) 



^((2/ -7/7 + 6)^2 +6(/7-2)^^ + 6ri)-(2n-;?) = (12) 
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\'S,Q{2n-p){2n-p-\)- 



(20p^ -\56 p'^ + All p^ -All p + \ 80) S"" + 
120(/ -6/7^ +ll/7-6)(J^^+ 

60(5/- 20/7 + 18)jVi + 
36Q{p-2)dyi'^ +\'^0y^ 



(13) 



Eliminating 5 from (12) and (13) gives two polynomial equations in ^ 

. Some lengthy algebra shows these polynomial equations yield the same quartic 
equation in which simplified gives: 



+24;7(;7-2«)^ 



210 + p 



4(;7-2«f (2;7-3)(30 + ;7(-50 + (37-9;7);7 + «(8;7-19)))^ 

-380 + «(l91 + /7(-303 + (l55-26;7);7)) 
+77(1 17 + /7(l39 + 2/7(9/7 -50))) 

^-(27 + 4/7(2/7-7))(-10 + /7(9/7-16)) 
« (-669 + 2/7 (596 + p{62p-335))) 

+I2p\p-2n){2p-l)(n + p{3p-n))n' +{l-2pf p\2p-3)n 



+A{2n-p)p 



-90 + p 



n 



=0 



(14) 



Although not repeated here because they are rather lengthy, it is immediately 
possible to solve (14) and hence write down formulae for the approximations to 
the first (xj) and last {Yp-\) zeros of ^{j^ , and the spacing parameter 5, 
directly in terms of any value of n and p. Using (7) then enables approximations 
of all of the zeros of Y(f ) to be obtained, thus enabling estimation of the spread 

and extent of i?„ / and the separation A?^ =(7^- -7;_2)« Z7q /2 of each of its 
zeros. 

Approximations are excellent for values of n and p that can be checked by direct 
numerical solution of the relevant polynomial. For example, even when the 
linearity condition, p«n and n»l is only marginally met, as when n-3Q 
and /7 = 10 , the approximation formulae zero set is {22.0, 27.6, 33.7, 40.5, 47.9, 
56.0, 64.6, 73.9, 83.8} while the direct numerical solution zero set of ^{y) is 
{22.4, 28.3, 34.3, 40.5, 47.3, 54.8, 63.2, 73.2, 85.9}, that is accuracies of -5% 
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or less. However when the linearity condition p«n and n»\ is well 

satisfied, the degree of accuracy of zero position approximation is remarkable. 
This can be seen in the following. For small sufficiently small values of p, the 
eigenfunctions may be exactly calculated analytically even when n»\, since 
the radial eigenflinction component contains only p terms. These analytically 
obtained values may then be compared with those obtained using the empirical 

relationship (7). With « = 1.0x10"''^ and p — ^, the present square root 

approximation treatment here gives the positions of the zeros of i?„ / (which are 

located at a radius r of -3.6x10^^ m) to better than 1 part in lO'** when 
compared with the numerical calculation of the zeros using the exact formula for 
i?„ „_4 . This leads to the accuracy of the actual spacing between the zeros (in 

this case ~ 44 km) being better than 4% even at these vast distances. 

3.3. Significance of the zeros of the radial eigenfunctions for overlap integral 
values 

When the initial and final wavefunctions do overlap, there are other reasons why 
the overlap integrals may still remain negligibly small in the case of large n. In 
dipole transitional decays originating from a 'deep' state (p » 1) and ending on 
a low p state, such as D to A of figure 1 where there is only unit change in A/ , 
the two states usually still have extensive overlap despite a large change in n 
because the much larger value of p in the initial state means that it can be spread 
over a large radial extent. One clear reason for a low value of in this case is 
that the normalisation condition on the upper eigenfunction requires that its 
average amplitude is very small because the function is highly spread. Another 
reason concerns the relative spatial oscillation frequencies (RSOF's) of the two 
states, particularly that of the high-p, upper state, which we now discuss. 

The preceding section enables us to determine the spacing of the zeros in the 
upper level and compare them with the width of the /» = 1 level. Considering 
again a central mass of 10^^ kg, the use of equation (14) shows that for 

A« = 10 (a proton quantum number change necessary to result in emission of a 
0.3 MHz photon) the spacing of those zeros in the higher-« eigenfunction that lie 

under the profile of the initial eigenstate is ~ 2 x 10~^ m, and constant (due to the 
particular value of the spacing parameter in these cases) over the whole of the 
radial eigenfunction at A in figure 1 to more than 1 part in 10^. This means that 

there are over -4x10^' uniform oscillations of the initial wavefunction under 
the relatively symmetric Pf - ^ final state function profile. The cancellation 

effect of this oscillatory behaviour on the radial dipole integral can be seen by 
noting the approximate value of the integral given by the formula 
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1 



^ 2 ^ 



23 



4 f e 1 2 

2 



4/ 



Ij^-XOr^''^^ bQ when Pf = l, «, «10'\ 



,34 



A« = /7,»10-^ [1]. 



4. Decay rates 

We now examine numerical values of, or upper limits to, the dipole decay rates 
of some of the gravitational eigenstates located at various positions on figure 1 . 
These numerical values are obtained either by direct integration of the relevant 
overlap integrals (small p), application of the various approximations developed 
in this series of papers or, in the most difficult cases, orders of magnitude 
estimated from the degree (or lack thereof) of overlap of the component 
eigenfunctions and their relative spatial oscillation frequencies. 

4. 1. Decay times for states along the p = \ diagonal ('Aw = 1^ 

The transitions involving only states along the j9 = 1 diagonal with high-« values 
have directly integrable radial overlap integrals, and potentially have the longest 
radiative lifetimes since only one dipole decay channel is available, the adjacent 
A« = l, = 1 state (B to A in figure 1). Table 2 gives dipole decay times for 
some of these An = 1 transitions for various central masses and relevant, suitably 
sized quantum numbers. The values are 'representative only' of the various 
transitions, in the sense that a nominal value of 0.5 has been taken for the 
angular component of the dipole matrix element. (Individual transition integrals 
between the various substates will vary considerably between values of this 
order and zero.) Their purpose is to give an indication of expected transition 
rates for various conditions and average upper limits on their values. 



17 



Gravitational eigenstates in weak gravity II: 



Table 2. Binding energies, positions, frequencies and decay times for high-«, 
p=\ eigenstates 



Central 
mass 
(kg) 


Eigenstate 
particle 
mass (kg) 


Quantum 
nuiiiDcr, n 


Energy 


Radius of 
eigenstate 
(m) 


Emission 
frequency 
(Hz) 


Dipole 
decay time 
(s) 


X 1 0^'* 


1 7x10~^^ 


10^^ 


91 


1 0x10^ 


~1x10~^ 


1 5x10^^ 


U.U ^ lU 


1 x10~^^ 


5 5 X 1 0^'* 

J.J ^ iU 


11x1 0~^ 


1 0x10^ 


~ 1 X 1 0~^ 


1 SxlO^^ 


9 0x1 0^^ 


1 7x10"^' 




6 9x1 0~^ 


1 0x10^^ 


~ 6x10"^^ 


~8x10-^^ 


2.0x10^'^ 


9.1x10"^' 


3.2x10^*^ 


3.7x10"^ 


1.0x10^^ 


-6x10 


-8x10^^ 


2.0x10^0 


1.7x10"" 


1.0x10^^ 


2.3x10^ 


3.0x10^ 


-1x10"^ 


-1x10^^ 


2.0x10^" 


9.1x10"^^ 


5.5x10^^ 


1.2x10"^ 


3.0x10'' 


-1x10"^ 


-Ixio'^ 


2.0 xlO"*^ 


1.7xl0~" 


1.0x10^"* 


2.3x10^ 


3.0x10^^ 


~lxlO~^^ 


-1x10^^ 


2.0x10"^^ 


9.1x10"^^ 


5.4x10^" 


1.3x10"^ 


3.0x10^^ 


~lxlO~^' 


-1x10^^ 



4.2. Other lm=l transition rates (transitions involving constant p) 

States not on the leftmost diagonal involve multiple allowed ( A/ = ±l) decay 
channels and we first consider the primary A« = 1 fransition decay channel of 
these states (that is transitions along constant p diagonals). Transition rates for 
An = 1 transitions along the second leftmost diagonal {p — 2) of figure 1 are 
relatively easily obtained by integration of the wavefunctions and yield similar 
values to those along p = \. To examine fransitions along higher p diagonals 

we need to specialise the general formula for the radial integral part /^j of n,j 

[!]./« (for / = — 1) may be re-expressed in terms of the final primary quantum 
numbers Uj- and Pf and the size of the quantum level jump A« as 



18 



Gravitational eigenstates in weak gravity II: ... 



r^2nf-lpf+J, I 2 A \' 

1 ' Irij +njAn\ 

Id = ~; 



nj-pj+2 



(/Jy +A«-2)! 



((2«y. -;?y)!(2«y. -;?y+A« + 1)1)2 X (15) 
^ {Irif -2pf +kf +2){2nf -2pf +kf +'}>){2nf -2pf +kf +A) ^ 



Pf-i 



kfKPf-kf-iy. 



(-2nf-2An) ' 



2Fi{a,b;c;z) 



where 2^1 ('^'^;<^;^) is the hypergeometric function with 



a = 2nf-2pj-+kj-+5, b = 2- p^ -An , c = 2nj- -2pj- +A and z = 



2«y + A« 



For the speciahsed case where A« = 1 and pj->2 , this may be written as 



_ (\ + 2nf)b, 
in — ~, r 



"~ 8(„,.l) 



pnf-pf^\)(2nf-pf^2)^ 



Pf 



n [2nf -2pf +1 + 2^ 



(16) 



!=2 



Z 



(2«y -2pf+kf + 2)(2nf -2pf+kf+ 3)(2«y - 2/7y + A:y + 4) 

(-2«.-2)*^ ^ ^ 

^ J— 2FAa,b;c;z) 

(2«/+l)'^ 



Provided />(^=/)y )««y and «/ »1 we expect that the decay times for all 

A« = l transitions along the lines of constant p = 2,3, 4... will be similar since 

the overlap is essentially 100% and their radial extent is limited, even though 
each successive pair of states has a different shape to the p = l pair (i.e. we 
expect that the A to B and E to F transitions of figure 1 will have similar decay 
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rates ). The hypergeometric function scales as order n ^ ^' , and it can be shown 
that this leads to + 1) / ~ / («y +1) . Thus for n^-X^'^, Ijf( 

Pj^ - I, Ap = 0, An - I) differs from for other p-pf^\ diagonals by less 

than 1 part in 1 0^"* per unit /(-diagonal change provided p « n . Hence 
corresponding transition rates for constant-p dipole transitions for deeper states 
such as E to F of figure 1 are essentially the same as for an A to B transition, 
provided p«nj- . It remains therefore to examine the transition rates of the 

other decay channels contributing towards the net decay rate of any given low-p 
state. These involve transitions where A« = 2,3,4... . 

4.3. Transitions where An = 2,3,4... 



We take equation (15) and modify it for large nj using Stirling's formula to 



avoid terms in nf\: 



5 

2' ^«/(«/+«/A«)^ (2«,.-;.^+l)^^^ h, 



\nf-Pf+2 



exp((/7^ -3)/2)(2«^ +AnY^ {rif -pf+lj' ''^ 

{Pf-l}. 



i^Pf+An-iy 



J. 

2 I &n+pj - 
X 

i=\ 



fAn+pj-2 

n (2nf-2pf+i + 3) 



V 



(17) 



kj-=0 



(2nf -2pf+kf + 2){2nf -2pf+kf + 3)(2«^ -2pf+kf+A) 

{-2nf-2Anff 
1 — 2F\[a,b;c;z) 



{2nf+Anr 

This equation includes terms of order n"^ which cannot be calculated when 
«y is large. Provided however that nj- » An and nf»pf, equation (17) may 
be put into a suitable approximate form, using the fact that 



lim 



, An 
1 + — 

V "/ 



2nf 



1 + - 



An 
2nf 



-lYf 



V 



which gives 
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— 



2Pf 



ipf-lMlnf-pf+lY"^ ^^Pf-2 



^2 



(2«y -2pj +kf + 2)(2«y -2pf + kf+ 3)(2«y -2pj + kf + A) 



kf\{pf-kf-\)\ 



{-2nf-2Anff 
{2nj-+Anf' 



(18) 



or equivalent alternate expressions involving various combinations of the other 
related parameters such as «, , p, etc. 



Equations like (18) not only enable easier approximate calculations of An = l, 
Ap-0 transitions than does equation (16), but also transitions where A« > 1 and 

ones which do not necessarily end on the Pf - ^ diagonal, provided that p f and 

A« remain small compared to . (The upper limit on p ^ for being 

calculable is determined by the number of terms one can tolerate in the 
summation and number of significant figures that are needed to avoid 

summation truncation errors.) The accuracy of the approximation can be 
estimated by the error associated with the above limit and that for Stirling's 

approximation, which were used in its derivation. (For the error in 

equation (18) is less than 1 part in 3000 when A« and /?j < 10 .) 



4.4. Decay rates for states where p^\ 



For the deeper {p>\) states of figure 1 with many decay channels, we designate 
the An-\ transition as the 'primary' decay channel. It might be expected that 
deep states could have so many channels available that the net decay rate could 
be much higher than that due to the primary An = 1 transition alone, especially 
given that these channels correspond to larger values of An which result in 
proportionally larger values of in a transition rate formula which involves a 
cubic dependence on <y. For a given deep state of fixed level n = lying on a 
diagonal p = p^ we wish to be able to calculate transition rates for each 
transition of the ensemble of A/ = 1 channels available, starting with An = l : 
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(nj,Pi)^ (Hf -Hj - l,pj- - pj) and continuing to An^p^: {n^,pj)^ 
(rif ^rif- Pj,pj- = \) (and an equivalent set for Al--\). For and values 

that might be significant in an astrophysical context (i.e. >10^^), Ij^ obtained 
using equation (18) is realistically calculable for p values up to several 
thousand). Figure 6 shows decay rates of all the allowed (A/ = l) channels for 
each of 9 initial states («, = 10^°,;?, = 25, 50, 100), ( = 10^^;7, = 25, 50, 100) 

and («, = 10^^,/?; = 25,50,100 ). The dramatic decrease in transition rate 1/ r 
demonstrates that the An — I is indeed the predominant decay channel of any 
high-H low-j3 state. 




An 

Figure 6. Decay rates (l/r) of all available channels of the initial low-p 
states Pi = 100 (short dashes), j3, = 50 (long dashes) and = 25 (solid curve) 
for an initial quantum level (a) = 1.0 x 10^"*; (b) n, = 1.0 x 10^°; and (c) n,; 
= 1.0 X 10^^; plotted as log(l/ r) versus the quantum jump An ; 
Parameter values: central mass = 2.0 x 10"*^ kg; particle mass = 1.67 x 10"^^ 
kg; charge = 1.6 X 10"'" C; 
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The decay rate of any particular channel is determined predominantly by the 
value of /jj as calculated by equation (18). Its use is limited to low p values, 
chiefly because of the number of significant figures required to avoid truncation 
errors, and the fact that the summation component requires ~p terms. However 
provided that n^ » 1 , estimates of channel decay rates may also be determined 
by considering an equivalent transition having the same Ap but one that ends on 
Pf=\ {c.f. equation (16) above). The value of /^j for these =1 transitions 

1 ft , , A -3 

2 



may be obtained from = ^q«^ 

» p^ » 1 , the logarithm of which is calculable for all Pj and 



2 

, valid when 



Figure 7 shows the variation of log(/g ) as a function of Ap= -l=(An-l) 

for dipole fransitions ending on =l) with «j remaining constant, as a 

function of the difference in the p value between the initial and final states. As 
expected, the values of Ij^ for low An are similar to their counterparts that were 
used to calculate the transition decay rates shown in figure 6. For the large 
values of (or «j ) considered here it can be shown from equation (16) that the 

ratio Ijf(p-2)/ Ij^(p-V) is ^J2n- . A similar trend continues: while p«n 
large, the ratio (/? + !)/ Ij^ (p) , is given approximately by (using equation 

(18)) i^~(4J72^){{p+i)/pY-' ~b,4^ I lUj and demonstrates that 

decreases dramatically (because of the large values of n^ and n^ ) for each 
successive unit increase in p . This results in the successive values for the 
overlap integral being also reduced by successive factors of ~ ^p I In^ for 
successive channels originating from the fixed upper state (because provided p 
is small, we can relate each decay channel from the given initial state to an 
equivalent, similarly-lifetimed fransition that ends on a p = \). 
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Ap(=A«-l) 

Figure 7. Variation of the radial component of the dipole matrix element 
plotted as a function of \o%{Ir) versus the difference in the p value between 
an initial {p^) and final state {pf=\).:nf = (solid curve); nf = 5 x 10"'^ 

(short dashes); nf= 10^"^ (long dashes). 

Parameter values: central mass = 2.0 x 10"^^ kg; particle mass = 1.67 x 10"^^ 
kg; charge = 1.6 X 10"''C; 



The trend Ij^ ~{^^p/ 2n- j((j9 + l)/ pY ' does not continue once p becomes a 
substantial iraction of Uj . Using equation (18) a better but more complex 



approximation than - b^n^ 



A-3 



e 
v2y 



for transitions ending on 



the p f - \ , and without the restriction n^ » p^ » 1 (although » 1 still a 
requirement) can be obtained as 

I 4>ly+2p,-+l 

y(2«^+/7,+i) 4 



_2"' Z?o(a +«/) ' 



nf+2 



7-2j 



(2nf+p,) 



Irif+pi+l 



(19) 
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where = ni + p , and from this equation it may shown that 
Ij^{p + \)l Ij^{p)^l as /j^ooso that the ln(/^) vs p plot of this fimction 
approaches a horizontal asymptote. 

Nevertheless, provided Pi (= A« - 1) remains small compared to n^ , the above 
discussion implies that the state-to-state fransition rate (1/r) for successively 
lower channels (i.e. as An increases) is <x:ay'Il^Jr^ oc A«''(A«/ 2«; )^ and also 

decreases rapidly. For example in the case of «, =10^\for an upper, moderately 
deep state (i.e. one whose initial state p value /?, is still much less than its n 
value «, ) the value of (and lijf ) for the Ap =1 fransition channel will be 

-10"'^ times that for the corresponding Ap = fransition channel, but ay" will 
increase only by a factor of 8, and the fransition rate for the Ap = 1 fransition 
will therefore be ~10^^ times smaller than for the 2^ = fransition. 

The net fransition rate for such a deep state («,,/, = - p) with » p»\ is 
then the sum of all the individual state-to-state fransition rates, 

p 

Z 4«i,«,-P).(«i-A«,«,-p-i) where 4„.,„.-p),(„,-A«,«,.-p-i) is the state-to-state 

A«=l 

fransition rate for the each of the respective decay channels. When is large the 
only significant term in this summation is the primary diagonal one 
An = 1, Af> = . Thus provided p « and » 1 , the lifetimes of any 
state whose parameter values are {n^,l^^n^- p) where n^»p>l, are 
approximately the same as the single state-to-state transition time of the 
primary diagonal channel for that state decay, that is the An = 1, Ap = 
fransition. This then is itself very similar in lifetime to that of the 
equivalent An = l, Ap = 0, /> = 1 state given in Table 2. 

If the condition on p«n, is not satisfied fransition decay rates increase 

dramatically with increasing An = l because of the factor co , and there is 
opportunity for a states to rapidly decay given the number of channels available. 
For example from figure 8 it can be seen that the value of (and hence the 

decay rate) for a given constant An = - «y increases to a non-trivial value 

once the initial p value is a significant fraction of n. The effect is sfrongly 



25 



Gravitational eigenstates in weak gravity II: 



dependent on (or equivalently n^) and the larger the greater the relative 
increase in on moving horizontally across the n-p state diagram. 
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high-/ 
transitions 



low-/ 
transitions 




p = n,-\ 



n = n 



/™=,=0 



8.0 



4.0 




Pf 



- 0: maximum angular momentum; 
= n^: zero angular momentum 



Figure 8. Variation of the radial dipole integral with the horizontal 

position of the transition on the n-p state diagram. An , A/ and n j- are 

constant for each curve, small dashes: «, = 1000, «/= 990; large dashes: 
«, = 750, «/= 740; solid: «, = 500, nf= 490. 
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The trend displayed here, that Ij^ has a significantly high value when p is a 

substantial fraction of n, means that decay rates of these low angular momentum 
states have potential to be very large. In contrast to the states where /> is a 
negligible fraction of n and multichannel decay is insignificant, this result 
demonstrates the critical importance of multichannel decay in low-/ cases (and 
also in low-«,- cases). Furthermore it enables us to understand the differences in 
interaction properties of the specialised, low-interaction, high-«, low-p states and 
those of the highly mixed states of more localised particles that have extremely 
broad eigenspectra, and how the Compton cross-section, for example, can be 
different for a set of such low-/?, high-« eigenstates compared to that obtained 
from the Klein-Nishina equation, from an ionised gas cloud or measured in the 
lab. 



4.5. Decay rates for states where p»\ 

When both initial and final states have p» \ then no approximations discussed 

so far can be used to estimate . It may still be possible however to gain 

information about the value of using knowledge of the position, extent and 

spatial oscillation frequencies of the state functions, found by solving equation 
(14). This could be particularly useful in particle interactions where one is 
concerned to know the range of states into which any given eigenstate might be 
scattered by interactions with a traversing particle, or the degree to which that 
traversing particle might itself be scattered. 

There is an interesting empirical relationship between the value of and the 
relative spatial oscillation frequencies (RSOF) of the two states involved in the 
overlap. Figure 9 shows that for transitions where one of the states has p = \ the 

value of log(-log(/^)) is linearly related to the log(RSOF) of the two states. 

This is a direct consequence of similar relationships that exist between p and /^j , 
and p and the spatial extent and spatial oscillation frequency of the wavefiinction 

shown in the insets. Thus given the p values of two states and the types of 
relationships shown in figure 9, it should be possible to predict their relative 
extents and relative spatial oscillation frequencies and from these obtain 
estimates of /^j for these states when they overlap. The significant thing about 
this approach to estimating values of I^ is that it has a much wider scope. 
Although it is a much more approximate method of determining upper limit 
estimates of overlap integrals, it has the advantage of depending only on the 
properties of the states themselves and not on the type of interaction potential. 
Thus it may be used as "first guess" or estimate to determine whether a decay or 
transition rate will be insignificant for other types of interactions as well, for 
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example non-dipole decays and particle interactions. If the spatial extents and 
oscillation frequencies of two states differ by large amounts then the probability 
of any interaction inducing transitions from one of these states to the other will 
be small even if the states overlap in space. 




15 

log(-log(/jj)) 

Figure 9. The value of log(-log(/;j)) versus log(RSOF) for 

transitions where one state has p = l and the other p = 4, 10, 10^, 10^ 
... 10^' (see text). Insets show the variation of the radial spatial extent 
(radial width) and RSOF with p. Parameter values: central mass = 2.0 
X lO'*^ kg; particle mass = 1.67 x 10"^^ kg; 



4. 6. Invisibility and the weakly interacting nature of the low-p eigenstates 

The extreme longevity of the low-p, high-« eigenstates when in large cenfral 
mass potentials has consequences beyond that of negligible electromagnetic 
radiation emission, for example their inability to effectively undergo 
gravitational collapse. More importantly however they will be weakly 
interacting with respect to a large portion of the electromagnetic spectrum. They 
will therefore be essentially invisible as may be shown from the following 
examples. 

We consider an incident flux of electromagnetic radiation impinging on an 
assembly of eigenstates and ask firstly, what transition rates does this radiation 
produce among the eigenstates, and secondly, what effects do the eigenstates 
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produce on the traversing beam of radiation? Ttie relevant processes for bound 
eigenstates are those of stimulated emission and absorption. (Classical 
electromagnetic scattering processes from 'free-particles' (Raleigh, Compton, 
etc.) may also be analysed using the summation of interactions over a vast 
eigenstate ensemble but the significant interaction channels for these processes 
do not involve channels that incorporate the low-/? eigenstates.) The probability 
per unit time /Jy of stimulated emission or absorption of radiation from a state / 

to a state / is given by [5] 

where p{cOif) is the radiation energy density per unit angular frequency and 
I(o}ff) is the corresponding beam intensity per unit angular frequency. 

Technically, for degenerate states, (20) must be summed over the available final 
states and averaged over the initial ones but since the significant transitions are 
dipole, the relevant selection rules mean that we may consider the effective 
degeneracies of the initial and final levels to be approximately the same and take 
(20) to be a reasonable approximation to the net eigenstate fransition rate in the 
following examples. 

As a first example, we consider photon absorption/stimulated emission from 
high-«, p = \ eigenstates. n,y is largest for A« = ±l (see equation (27) of [1]). 

The most pervasive radiation extending over the entire volume of the eigenstates 
is that of the CMB and at the resonant frequency at present day CMB densities, 

yO^ta^y ) ~ 7x10"^^ Jm^^srad ' . Taking a typical proton-occupied eigenstate 

with «~5xl0^^ in a gravitational potential of central mass -lO'*^ kg gives 
i^~4xl0^ s We thus expect that any typical = ! eigenstate to be 
subjected to random promotion/demotion state fluctuations of this magnitude. 
Over a cosmic time of 10^' s, any given state might be expected to have 

dispersed over \/4xl0^^ -lO'** w-levels, clearly negligible compared to the 
value of n itself Even allowing for an increased intensity of >o(<y,y ) over cosmic 

history by using /?(ft>^) ~ 7x10^* Jm"^srad ' (at z = 10000) gives a «-level 

dispersion of ~ 10^^ over the cosmic history. Similar values will be obtained for 
other low-/?, ^ = , fransitions. 

There are of course many promotion channels (A« >1) available for excitation 
of p = \ (or other low-p) eigenstates. However, as demonstrated by equation 
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(27) of [1] and equation (18) of this present paper, the value of the radial part of 

Tljf decreases so rapidly with increasing A/t , that the contribution to the 

excitation of a low-p eigenstate by stimulated absorption is dominated by the 
Ap = channel. 

The second example involves calculation of the numbers of photons that would 
be absorbed or scattered in crossing a typical halo. We consider a beam of 

photons of number density «^ m whose frequencies are distributed uniformly 

within a small angular frequency range Aty . The beam intensity per unit angular 
frequency I{a>) is /(<y) - n^cticol Aco , enabling (20) to be rewritten as 

P..= (21) 

where P^^ represents the probability per unit time for the promotion of a single 
eigenstate. We assume that n^ is uniform throughout a typical halo. The number 
of excitations per second within the halo is then given by P^j times the number 
of available eigenstates A^^ that have transitions whose frequencies lie within 
the angular spread Aa> . For the purposes of illustration we assimie that the halo 
consists of an eigenstructure filled from the longest lived p = l states to some 

maximum p level . To accommodate a halo mass of -10"*^ kg with radius 
-l.SxlO^^m requires totally filling all levels from p = l to p = \0. 

We now consider the transition An-\,Ap -0 , originating from the state 

n = 5xlO^^ p = \. For a given value of n and p, there are 2/ + l~10^'* z- 

projection substates available for promotion with the same central frequency. It 
may also be shown that Ap = fransitions originating on the states 

« = 5xlO^^, /7 = 2,...10 have profiles that lie essentially within the natural 

linewidth of the « = 5xlO^^^ P = ^, An = l,Ap = fransition, and that so do 

states whose «-values lie within ~ 1 x 10^^ of « = 5 x 10". This means that if we 

take the spread Aco of the riy m~^ photons to be within the natural linewidth of 

the transition, there will be 10"^"^ x 10x10^^ = 10^^ states available for promotion 
within the angular frequency spread of the photon beam, so that ~10^^. 
Since every time a fransition is excited by absorption a photon is removed from 
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the field we may then, using (21) write down a rate dN^ I dt at which photons 
are removed from the photon beam across the whole halo as 

dN^ Att^IIJ coN, 

— r_ = i — r_ i (22) 

dt SsohAcD 

or by writing = n^V where Fis the halo volume, obtain the expression 
dn^ = -(1 / t) Hydt , where t is the expected lifetime for a photon of frequency 

CO to survive without absorption, given by 

r= y^^" (23) 



For a transition An = l,Ap = 0, originating from « = 5x10 p = \,T~ 10 s 

so that photons of frequency ~ 10~^^ Hz would be readily absorbed as would be 
expected from the first example. For An = 2,Ap = l, originating from 

« = 5 X 10^^ ^ /> = 1 , the absorption frequency is twice as great but Tiff is much 
reduced so that r ~ 1 0~^^ s . For M = 3,Ap = 2 the hfetime is of the order of the 

galactic halo fransit time ( r ~ 10^^ s ) and for A« > 4, Ap = 2 the photon lifetimes 
become much greater than the halo transit time. Thus for any photon frequency 
realistically observable within the electromagnetic spectrum ( co>l kHz implies 

1 8 

A«>10 ), p = l eigenstates do not absorb radiation. Provided the p-value of 

the initial state is small, and certainly for />< 1 0^ , similar results would apply, 

leading to the invisibility of all but the highest /7-valued sates. Therefore 
provided that an eigenstructure consisted predominantly of particles occupying 
\ow-p eigenstates, it will be transparent to virtually all regions of the 
electromagnetic spectrum including the optical region. 

The weak response of the low-p eigenstates to stimulated emission and 
absorption comes about directly as a consequence of the rapid decrease in the 
value of with increasing RSOF between initial and final eigenstates, 
discussed in the previous section. Similar cancellation effects within the overlap 
integrals can be seen in the angular domain as well and whenever Ap, An, Am 
or A/ » 1 there can be significant differences in spatial oscillation frequencies 
between initial and final states resulting in weak interaction often irrespective of 
the Hamiltonian involved. Given an initial high-«, low-p state on figure 1, one 
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can in general one can make a distinction between several types of transitions 
depending on the location of the final state on figure 1. In the case of photon 
absorption discussed above for example, the final state lies almost vertically 
above the initial state (since A/ = ±l) and there is always a large change in p 
giving rise to large RSOFs and infinitesimal values for . We can also consider 

interactions giving rise to large horizontal state changes. This corresponds elastic 
scattering: large changes to the momentum of both the eigenstate and the 
perturbing entity occur without energy transfer. Necessarily large changes in p 
are involved with correspondingly low overlap integrals and hence limited 
opportunity for interaction unless the selection rules are such as to permit a large 
range of available channels. It is clear that when processes such as elastic and 
Compton scattering are described in the present basis set of eigenfunctions, then 
these processes must be seen to take place largely through mixtures of the low-/, 
high-/7 channels that enforce the traditional localization the particles and give 
them opportunity to interact in the classical sense. 

To avoid changes in p and ensure non-negligible overlap integrals, interactions 
need to involve initial and final states that lie on relatively closely spaced 
diagonals. For low-/? states this presents a problem for interaction. The fact that 
the radial eigenfunctions decrease rapidly outside their range, coupled with their 
relatively small radial spreads and positions that vary significantly with n, 
enforces a limit on the extent to which momentum and energy can be 
transferred. If the initial and final states differ too much in energy and/or 
momentum then they will cease to overlap and interaction integral will be zero, 
irrespective of the potential involved. This is major reason for the weakly 
interacting nature of the low-p states. 

Equation (14) is useful in determining when two given states will cease to 
overlap and in defining the menaing of 'low-/?'. For example given two deeper-p 
states (but still with the condition that p«n) we may ask a similar question 
that posed earlier for transitions between two p = \ states, that is: how far away 
must such states be on the n-p diagram to ensure that 'no' transitions occur 
between them irrespective of the interaction potential involved? Essentially this 
involves comparing the spatial position of the lowest zero of the radial function 
of the upper level with the spatial position of the largest zero of the radial 
function of the lower level. If there is a substantial positive difference between 
these two values then we can be reasonably confident that the radial functions do 
not overlap because the mathematical behaviour of the radial eigenfunction is to 
decrease rapidly for values outside its radial range. This work is still in progress 
and will be reported on at a later date. We do expect however that cross sections 
for interactions between particles and eigenstates will be larger than those with 
photons, because of the less restrictive selection rules for interactions involving 
particles. (That is low-p eigenstates may not be as "invisible" for particles as 
they are for photons.) 
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We emphasise one final point about ttie measurement of cross sections. The 
classical elastic or inelastic photon scattering processes (for example Compton 
and Rayleigh scattering) have measured cross sections determined from 
experiments carried out using particle assemblies in vessels that contain 
macroscopic numbers of particles localized within the vessel. These 'localized' 
particles necessarily have eigenspectra that, by virtue of this localization, are 
composed of a vast but reproducible summation of the much more highly 
interacting low-/, high-p states, the makeup of which is determined by 
macroscopic parameters such as temperature, number density etc. Individual 
eigenstates or specialised arrays of eigenstates with unique features will have the 
capacity to exhibit significantly different properties from say, a gas in thermal 
equilibrium (in the same way for example that the scattering cross section for a 
atomically bound electron differs from that of a free elecfron). Any measurement 
of a particular cross section should necessarily be accompanied by reference to 
the relevant time-averaged eigenspectra of the particle ensemble under which the 
measurement was made. This is done automatically for example when we 
recognise the dependency of cross sectional data on external parameters such as 
temperature etc. that determine the eigenspecfral mix. 

5. Conclusion 

This paper has developed further, more approximate procedures for studying the 
high-n, high-/, low-p gravitational eigenstates. This has included investigation of 
the trends in functional behaviour of both the angular and radial components of 
high quantum number eigenfunctions and development of a formula for 
calculating the width, position and spatial oscillation frequencies of the zeros of 
the radial component. These trends have enabled many potential dipole decay 
channels to be eliminated as effective channels and hence facilitated estimation 
of the net dipole decay rates of many 'deep', high-« states even when those states 
may have many channels available for decay. 

The first point of interest is the fact that the spatial oscillation frequencies of the 
eigenfrinctions increase in proportion to the square root of the /j-value, which is 
the dominant physical explanation as to why the radial overlap integral for high- 
p with low-p states decreases so rapidly with increasing p, irrespective of the 
interaction potential and the degree of state overlap. It means that not only do the 
low-/* eigenstates have extremely long decay times and gravitational stability, 
but also that they cannot easily be pushed into higher levels by stimulated 
transitions or (to a lesser extent) transferred to other locations with very different 
p values on the state diagram by other interactions such as particle collisions. 
This accounts for their invisibility and their generally weakly interacting nature. 
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An important consequence of the dramatic effects of spatial oscillation 
frequencies is that high-n transitions must necessarily take place along diagonal 
lines of approximately constant-p. For the high-«, \ow-p states, the only 
transitions possible are those whereby they are transferred to other low-p states. 
For these types of transitions however, the role of the distinct width, position and 
rapid cut-offs of the states becomes important. This restricts the amount of 
energy («) and momentum (/) that can be exchanged in any interaction involving 
a low-/? state and consequently also restricts the changes that take place in the 
perturbing entity and again contributes to the eigenstate's weakly interacting 
behaviour. On the other hand the low-/ states have large radial extents and 
interactions that approximately conserve p can take place with a much greater 
exchange of energy and momentum. This explains how the cross-sections for 
traditional processes like elastic and Compton scattering can be quite different 
for the low-p eigenstates. 

It is also interesting to note that the weakly interacting properties of the particles 
occupying the low-p eigenstates arise as a consequence of the quantum states in 
which they are placed and not from any internal weakly interacting properties of 
the particles themselves, a fact that is highlighted by the present work. Indeed 
because the weak interaction is a consequence of the properties of the 
wavefunctions themselves and their behaviour on the scales and sizes of the 
quantum parameters involved, the phenomena is not restricted to gravitational 
potentials and general macroscopic low-p eigenfunctions could in some sense 
provide us vwth an unusual structural form of matter that we are yet to become 
familiar with. 

The prediction of invisibility and weak interaction of the speciaUsed 
gravitational eigenstates is undoubtedly the most significant result to come out 
of the present work. If such gravitational eigenstates exist then they may be 
responsible for most of the dark matter in the universe, especially if, as it 
appears, a viable primordial-nucleosynthesis-consistent scenario for their 
production might be possible [2]. This would make them particularly 
attractive as an explanation for dark matter, since it means that its existence 
could be explained without the need to introduce any new particles and that its 
origin follows as an automatic prediction from fraditional quantum theory. 
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